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Abstract 

We study the antiferromagnetic g-state Potts model on the square lattice 
for q = 3 and q = 4, using the Wang-Swendsen-Kotecky Monte Carlo algorithm 
and a new finite-size-scaling extrapolation method. For q = 3 we obtain good 
control up to correlation length £ ~ 5000; the data are consistent with £(/?) = 
Ae 2 P /3(1 + a\e~^ + . . .) as fi —> oo. For q = 4 the model is disordered even at 
zero temperature. 



PACS number(s): 05.70.Jk, 64.60. Cn, 64.60. Fr 



The Potts model [1,2,3] plays an important role in the general theory of criti- 
cal phenomena, especially in two dimensions [4,5], and has applications to various 
condensed-matter systems [2]. Ferromagnetic Potts models have been extensively 
studied over the last two decades, and much is known about their phase diagrams 
[2,3] and critical exponents [5,6]. But for antiferromagnetic Potts models, many 
basic questions remain open: Is there a phase transition at finite temperature, and 
if so, of what order? What is the nature of the low-temperature phase? If there is 
a critical point, what are the critical exponents and the universality classes? 

In this Letter we report the results of a large-scale Monte Carlo study of the 
3-state and 4-state antiferromagnetic Potts models on the square lattice, using 
the Wang-Swendsen-Kotecky (WSK) [7,8,9] cluster algorithm. We use a powerful 
new finite-size-scaling (FSS) extrapolation method [10,11] to estimate the infinite- 
volume correlation length £ and staggered susceptibility Xstagg- Using lattices up 
to 1536 2 , we can attain an accuracy of a few percent on £ and Xstagg at correlation 
lengths £ as large as 5000. This allows us to conjecture the exact form of the critical 
behavior for the 3-state model. 

The g-state Potts model is defined by the reduced Hamiltonian 

n = -J£<W a , (i) 

(xy) 

where the sum runs over all nearest-neighbor pairs of lattice sites, and each spin 
takes values a x £ {l,2,...,g}. The antiferromagnetic case corresponds to J = 
-f3 < 0. 

Baxter [4,12] has determined the exact free energy (among other quantities) on 
two special curves in the (J, g)-plane: 

e J = 1±V? (2) 
e J = -l±^~q (3) 

Curve (2_|_) is known to correspond to the ferromagnetic critical point, and Baxter 
[12] conjectured that curve (3+) corresponds to the antiferromagnetic critical point. 
For q = 2 this gives the known exact value; for q = 3 it predicts a zero-temperature 
critical point (J c = — oo), in accordance with previous belief [13,14]; and for q > 3 
it predicts that the putative critical point lies in the unphysical region (e Jc < 0), so 
that the entire physical region — oo < J < lies in the disordered phase. 

We remark that the q = 3 model is exactly soluble at zero temperature in an 
arbitrary magnetic field [14,15,16]. We would not be surprised if the model were 
found to be exactly soluble also at nonzero temperature (in zero field), at least in 
the sense of determining the exact asymptotic behavior as f3 — > oo (see [17,18,19] 
for some tantalizing ideas in this direction). 

Nightingale and Schick [20], using a phenomenological-renormalization method 
based on infinite strips of width 2-8, confirmed the prediction of a zero-temperature 
critical point for q = 3, and claimed that the correlation length diverges as £ ~ 
exp(c/3~ 1 ' 3 ). Wang, Swendsen and Kotecky [7,8], using Monte Carlo, claimed to 
confirm this latter behavior. But this behavior seems a priori implausible to us: 
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the fundamental variable in the Potts model is t = e , so an ordinary power-law 
critical point £ ~ (t — t c )~ v with t c = would correspond to £ ~ e^. Moreover, we 
suspect that this model can be exactly solved, in which case v would most likely be a 
rational number. We are unable to imagine any mechanism leading to £ ~ exp(c/3 K ) 
with k ^ 1. 

We simulated the models for q = 3 and g = 4. We measured the energy E, the 
staggered susceptibility Xstagg, and the (finite- volume) second-moment correlation 
length 

l/2 

c I (Xstagg / ' F stagg) ~ 1 \ /.x 

\ 4 sin (vr/L) / 

where Xstagg = G(7r, 7r) and F stagg = G(tt + 2tt/L } tt) [here G is the Fourier transform 
of the spin-spin correlation function]. 

For q = 3 we ran onlxi lattices with L = 32, 64, 128, 256, 512, 1024, 1536 at 
153 different pairs (f3,L) in the range f3 < 6.0 (corresponding to £ < 20000). Each 
run was between 2 X 10 5 and 2.2 X 10 7 iterations of the WSK algorithm, and the 
total CPU time was about 2 years on an IBM RS-6000/370. The WSK algorithm 
appears to have no critical slowing-down: we found T int ^ M 2^ < 5 uniformly in f3 
and L [21]. We extrapolated to infinite volume by the method of [10,11], taking 
imin = 10 and L m i n = 64 and using a quartic fit for the FSS functions (see [11,22] 
for more details). We thereby obtained £ to an accuracy of about 1% (resp. 2%, 
3%, 5%) at ^ 1000 (resp. 2000, 5000, 10000). The errors on X stag g were roughly 
twice as big. 

Our data are in clear agreement with the prediction of a critical point at f3 = oo. 
The correlation length £ rises roughly like e 2 ^, and we initially thought that this was 
the exact asymptotic behavior. However, at /3 > 3.4 (£ > 75), £ begins to rise faster 
than this (Fig. la). We therefore guessed a multiplicative logarithmic correction, 
i.e. £ ~ e 2 ^ f3 p for some power p; see Fig. lb,c for p = 1/2 and p = 1, respectively. 
It is difficult to distinguish between 1/2 < p < 1 without additional information on 
the corrections to the leading asymptotic behavior. We do not know how to carry 
out a low-temperature expansion around the (critical) zero-temperature state; but 
it is reasonable to expect that there exists an expansion in powers of , which 
corresponds to a minimum energy cost of one unit for an "overturned" spin. We 
therefore expect 

C(/3) = Ae 2fj ^ [l + ai e~ fj + a 2 e~ 2fj + ...]. (5) 

If we accept this Ansatz, a value p 1 is clearly favored, with A 0.184 and 
a t w 11 (Fig. 2). 

The critical exponent v = 2 found here corresponds to an operator with scaling 
dimension X = 2 — 1/^ = 3/2, which is one of the possibilities proposed by Saleur 
[18, p. 248]. A logarithmic correction f3 p ~ (logt) p withp integer (particularly p = 1) 
can occur as a result of "resonance" between operators whose scaling dimensions 
are rationally related [23]. We hope that the numerical results presented here will 
serve as useful clues toward the exact solution of this model. 
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It is only fair to note that our data can also be fit reasonably well by £ ~ exp(c/3 K ) 
with n ~ 1.3 (Fig. 3). It is on theoretical grounds that we prefer (5). We doubt 
that these two Ansatze can be reliably distinguished by numerical means in the 
foreseeable future. 

The staggered susceptibility is consistent with the believed exact behavior [24] 
Xstagg ~ C 5 ^ 3 ? unmodified by any further powers of f3. The energy per site is consis- 
tent with 

E(/3) = b ie - fj + b 2 e- 2fj + ... (6) 

with b x w 0.22 (Fig. 4). 

For q = 4 the story is very brief: simulations on L = 32, 64 show that £ < 2 
uniformly as f3 — > oo. Clearly there is no critical point in the physical region. 
Physically, there is so much entropy that the correlations decay exponentially even 
at zero temperature. This can be proven rigorously to occur on the square lattice 
for q > 8 [25,22], and our simulations confirm Baxter's [12] prediction that it occurs 
in fact for q > 3. 

Details of this work will appear elsewhere [22]. 

We wish to thank Chris Henley, Paul Pearce, Andrea Pelissetto and Hubert 
Saleur for helpful comments. The authors' research was supported in part by the 
Conselho Nacional de Desenvolvimento Cientifico e Tecnologico-CNPq (S.J.F.), the 
Fundacao de Amparo a Pesquisa do Estado de Minas Gerais (S.J.F. and A.D.S.), 
and U.S. National Science Foundation grant DMS-9200719 (A.D.S.). 
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d=2, q = 3 AF Potts Correlation Length 
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Figure 1: Infinite- volume correlation length £ divided by e 2 ^ f3 p for (a) p 
p = 1/2, (c) p = 1. Error bars are one standard deviation. 
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d = 2, q=3 AF Potts Correlation Length 
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2: il{e 2l3 f3 v ) with p = 1, plotted versus e t Note the nearly linear behavior, 
good agreement with (5). 
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d = 2, q=3 AF Potts Correlation Length 
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Figure 3: Log-log plot of log£ versus f3. The indicated asymptotes are, from left to 
right, £ = exp(0.5469/3 1 - 7 ) and £ = exp(0.8465/3 1 - 37 ). 
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d = 2, q = 3 AF Potts Energy 
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Figure 4: Energy per site E divided by e ^, plotted versus e t Symbols indicate 
L = 32 (+), 64 (x), 128 (□), 256 (O), 512 (o), 1024 (*), 1536(*). Error bars are 
invisibly small. The uppermost points at each f3 represent the infinite- volume limit. 
Note the nearly linear behavior, in good agreement with (6). 
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